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Abstract Following recent work on the combination of elec- 
troweak and strong radiative corrections to single W -boson 
hadroproduction in the POWHEG BOX framework, we gen- 
eralize the above treatment to cover the neutral current Drell- 
Yan process. According to the POWHEG method, we com- 
bine both the next-to-leading order (NLO) electroweak and 
QED multiple photon corrections with the native NLO and 
Parton Shower QCD contributions. We show comparisons 
with the predictions of the electroweak generator HORACE, 
to validate the reliability and accuracy of the approach. We 
also present phenomenological results obtained with the new 
tool for physics studies at the LHC. 

Keywords Hadron colliders • Drell-Yan ■ QCD ■ Elec- 
troweak corrections 
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1 Introduction 

The majority of the searches and studies at the proton-proton 
(pp) collider LHC at CERN, including the analyses related 
to the Higgs boson, is based on an intensive use of Parton 
Shower (PS) generators matched with fixed-order pertur- 
bative calculations. The success of these general computa- 
tional frameworks, like MC@NLO [1] and POWHEG [2, 3], 
rests on their capability to provide reliable predictions for 
both inclusive cross sections and distributions of a large va- 
riety of signatures in the presence of arbitrary event selec- 
tion conditions. The increased theoretical accuracy with re- 
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spect to traditional PS generators has promoted the wide 
use of these frameworks on the experimental side. Gener- 
ally speaking, the processes available in such tools include 
finite-order contributions at the next-to-leading order (NLO) 
accuracy in QCD [4]. Only recently the implementation of 
electroweak (EW) contributions in matched PS generators 
has been addressed in the literature [5, 6] in the process of 
single W -boson hadroproduction, according to the matching 
approach given by the POWHEG method. High precision 
in the W hadroproduction process is required in particular 
by the need to reduce theoretical uncertainties in W mass 
measurements. The two implementations described in Refs. 
[5, 6] are both available as subprocesses in the repository 1 
of the general computer framework POWHEG BOX [7]. 

For many physics studies at the LHC, also the process 
of lepton pair production in hadronic collisions, known as 
neutral current Drell-Yan (NC DY), requires high theoretical 
accuracy. The state of the art of fixed-order calculations is 
encoded in the next-to-next-to-leading order (NNLO) QCD 
programs DYNNLO [8] and FEWZ [9, 10], in the EW codes 
HORACE [11, 12] and ZGRAD/ZGRAD2 [13, 14], while 
the SANC framework [15-18] and the RADY code [19] al- 
low to evaluate both QCD and EW NLO corrections. Re- 
cently, NNLO QCD corrections were combined with NLO 
EW contributions in FEWZ [20]. 

All these theoretical efforts were made because, given 
the high-precision measurement of the Z-boson mass at LEP, 
the NC process is very helpful, if not unavoidable, for de- 
tector calibration purposes at LHC and is a standard can- 
dle that can be used to constrain the Parton Distribution 
Functions (PDFs). Moreover, the transverse momentum dis- 
tribution of the Z boson can be accurately measured, and 
used to tune non-perturbative parameters in the generators. 



1 See http://powhegbox.mib.infn.it for an updated list of all available 
processes. 
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This measurement indirectly constrains the transverse mo- 
mentum distribution of the W boson, with an accuracy that 
strictly depends upon the accuracy of the production model. 
The NC DY process also allows to perform tests of the SM 
at the loop level and to measure EW parameters, such as the 
weak mixing angle, in a hadronic environment. In the high 
tail of the transverse momentum and invariant mass distribu- 
tions of the produced leptons, the NC DY is one of the main 
irreducible backgrounds to the searches for new particles at 
the LHC. 

These physics motivations require the simultaneous con- 
trol of all the relevant QCD and EW higher-order contribu- 
tions to the NC DY process. Following the recent implemen- 
tation of EW corrections to W production in the POWHEG 
BOX as described in Ref. [5], here we present the analo- 
gous treatment for the NC DY process. In particular, the phe- 
nomenological importance of combining QCD and EW cor- 
rections to dilepton hadroproduction is illustrated through 
the analysis of various numerical results. Let us note that our 
approach is complementary to the combination of NLO EW 
corrections to NC DY with fixed-order QCD at NLO and 
NNLO accuracy as recently realized in the SANC frame- 
work [17, 18] and in the FEWZ code [20], respectively 2 . 
However, unlike the above combinations, QCD and QED 
shower effects are taken into account in our realization. Thanks 
to it, it is therefore possible to obtain reliable predictions in 
the presence of combined QCD and EW effects also for ex- 
perimentally relevant observables, such as the Z and lepton 
transverse momentum, which are not realistically accessible 
to the tools including fixed-order corrections only. 

The paper is organized as follows. In Section 2 we de- 
scribe the modifications applied to the POWHEG BOX for 
the inclusion of the NLO EW and higher-order QED correc- 
tions in the NC DY channel. In Section 3 we present and dis- 
cuss several numerical results obtained with the new tool at 
LHC energies, both as a cross-check of the EW corrections 
at NLO accuracy and about the interplay of QCD, QED and 
weak corrections. Conclusions are given in Section 4. 



2 Details of the calculation 

POWHEG (POsitive Weight Hardest Emission Generator) 
is a method conceived for embedding NLO QCD compu- 
tations into PS simulations. Here we just review the basic 
ingredients and ideas of the method, paying particular atten- 
tion to the components that we generalized for the inclusion 
of the EW and QED corrections. For more details the reader 
is referred to the original literature [2, 3, 7]. 



2 Further studies on combining QCD and EW corrections to W/Z 
hadroproduction are described in Refs. [21-27] 



2.1 The POWHEG method 

In the POWHEG formalism, the generation of the hardest 
emission is performed first, using full NLO accuracy, and a 
shower Monte Carlo is used to generate subsequent radia- 
tion. Therefore, the building blocks of the method are those 
typical of a NLO calculation and of PS generators. 

At NLO accuracy in QCD, the necessary ingredients are 
given by the process-dependent virtual corrections Vqcd and 
real radiation matrix element(s) Rqcd- To ensure cancella- 
tion of the initial-state collinear singularities in hadronic col- 
lisions, further components of the calculation are two factor- 
ization counterterms, one for each of the incoming partons 
(©,0), which we denote, following the POWHEG notation, 
as G® C q (the collinear remnants). Soft and collinear diver- 
gences coming from real radiation are treated in POWHEG 
using a subtraction method (the FKS subtraction formal- 
ism [28]) which requires the introduction of a set of func- 
(a) 

tions Cq CD (the real counterterms), where the superscript a 
denotes a particular singular region of Rqcd- Correspond- 
ingly, the method requires that the real radiation contribution 
is separated in all those terms Rq CD that are only singular in 
a single region a. 

With all these ingredients at hand, the calculation of an 
inclusive cross section with NLO QCD accuracy proceeds in 
the POWHEG BOX according to a high degree of automa- 
tion. First the algorithm identifies all the singular regions 
and maps the real radiation matrix element over the singular 
configurations. Afterwards, it performs the subtraction pro- 
cedure and computes the collinear remnants. As shown in 
the following, this procedure can be effectively generalized 
to the inclusion of EW contributions, provided the appropri- 
ate ingredients and modifications are supplied. 

In the POWHEG method, the matching of the NLO com- 
putation with the PS is achieved in terms of a modified Su- 
dakov form factor, which contains the NLO real radiation 
matrix element and is equal to the product of the Sudakov 
form factors for each singular region. The generation of the 
event with the hardest radiation, which is a fundamental as- 
pect of the approach, is also handled automatically by the 
POWHEG BOX framework and can be adapted to cope with 
the radiation of colorless partons, as discussed in the follow- 
ing. 

2.2 Electroweak contributions: calculation and inclusion in 
the POWHEG BOX 

We treated the NC DY process according to the same recipe 
adopted in Ref. [5] for the implementation of the EW con- 
tributions to single W production in the POWHEG BOX. 

The most intricate and lengthy step is the calculation of 
the one-loop EW corrections Vew> since it involves the treat- 
ment of unstable particles in addition to the cancellation of 
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mass singularities. We accomplished this task realizing a li- 
brary for the computation of Vew, based upon the one-loop 
structure for the virtual photonic and weak corrections given 
in Ref. [19]. Unlike Ref. [19], where mass regularization is 
used to treat soft and collinear singularities, we employed 
a mixed scheme in order to comply with the use of dimen- 
sional regularization of POWHEG. Precisely, we used di- 
mensional regularization to treat the singularities associated 
to the QCD partons and the photon, but we kept the mass of 
the leptons finite, since this mass has a well-defined physical 
meaning, being the true regulator of the QED mass singular- 
ities. For the evaluation of the scalar integrals, we resorted in 
particular to various results of Ref. [29] (for the computation 
of the three-point functions) and Ref. [30] (for the computa- 
tion of the four-point functions). 

The singularities associated to the unstable nature of the 
weak bosons circulating in the loops and given by loga- 
rithms of the form log(s — My + is) were treated according 
to two different schemes, which both respect gauge invari- 
ance. In the first approach (factorization scheme) [19, 31], 
the logarithms that are singular at the resonance (and would 
be cured by a Dyson resummation of the self-energy contri- 
butions inside the loop diagrams) are treated with the substi- 
tution log (s — M v + ie ) log (s — My + ir v M v ),V = W,Z. 
In the second approach (complex mass scheme) [32, 33], the 
squared vector boson masses are taken as complex quantities 
jX wz = M^z ~ itw,zMw,z in the LO and NLO calculation, 
leading to complex couplings as well. Both schemes were 
implemented in the POWHEG BOX and the differences be- 
tween the two procedures for treating the resonance are at 
the per mille level or below it for arbitrary partonic centre of 
mass (cm.) energies s [19]. 







Fig. 1 The photon bremsstrahlung diagrams of the NC DY process. 

To complete the NLO EW calculation, we computed the 
real photon matrix element Rew relative to the bremsstra- 
hlung processes qq y,Z — ^— > l + l~y shown in Fig. 1. We 
included finite lepton mass effects in the calculation and re- 
sorted to the algorithm developed and described in detail in 
Ref. [5] for the treatment of the radiation phase space with 
massive fermions. The search for the singular regions was 



extended to handle the divergences associated to photon ra- 
diation. The original algorithm applied to the initial-state ra- 
diation (ISR) of QCD partons was generalized to treat the 
mechanisms of QED ISR and final-state radiation (FSR). As 
the large logarithms associated to photon ISR have to be 
reabsorbed in the PDFs, in analogy to QCD, the dominant 
QED contribution is due to FSR, which contains logarithms 
of the form a &m \og(s/m 2 ) for the typical non-fully inclusive 
and non-infrared safe conditions used in the event selection. 
As demonstrated in Ref. [5], the new phase space treatment 
including massive leptons is a key ingredient to take under 
control the above logarithmic contributions. To ensure can- 
cellation of all the mass singularities in terms of the new 
phase space, we also recalculated the QCD matrix element 
Rqcd including finite lepton mass contributions in the cal- 
culation of both the gluon emission qq —> y,Z l + l~g 
and gluon-induced gq/q — > l + l~q/q processes. 

To summarize, the implementation of EW corrections 
in the POWHEG BOX was realized through the following 
generalization of the native QCD elements 
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and improving the algorithm for the identification of the sin- 
gular regions associated to real photon radiation. 

For the cross section calculation, we implemented three 
input parameter schemes: the G^, the on-shell a em (Q) and 
the a em (Mz) scheme. In the first scheme, the primary in- 
put parameters are given by the muon decay constant and 
the weak boson masses; in the second one, is replaced 
by the fine structure constant at zero momentum transfer. 
In the latter scheme, the non-perturbative hadronic contri- 
bution to the running of a sm is treated in terms of effective 
quark masses reproducing the correct value of a em (Q 2 ) at 
high energies, i.e. Q 2 = M§. In the third scheme, the input 
parameters are oc em at the scale Q 2 = M%, Mz and Mw- The 
widths of the weak bosons are kept fixed everywhere. 

As a last step, we matched the NLO QCD and EW cor- 
rections with PS contributions. We considered both QCD 
and QED parton cascade. According to POWHEG, the gen- 
eration of the hardest-radiation event is performed in terms 
of the exact real radiation matrix element and of a modified 
Sudakov form factor including it. Once the configuration 
with the hardest (transverse momentum p±) emission has 
been generated, the subsequent radiation processes handled 
by the PS take place at lower p^, applying a veto technique. 
The above algorithm was generalized to include photon ra- 
diation in the POWHEG BOX. We made use of both R QCD 
and Rew and implemented a Sudakov form factor equal to 
the product of the form factors for ISR gluon and photon 
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radiation, as well as FSR photon radiation. In the native 
QCD construction, the method requires a lower cut-off on 
the transverse momentum, in order to avoid to reach un- 
physical values of the strong coupling constant and of the 
PDFs. Also this requirement was generalized, following the 
procedure applied to W production in Ref. [5]. Precisely, the 
lower p± cut-off was set equal to a typical hadronic scale for 
gluon or photon radiation from quarks, while it was taken as 
the mass of the lepton for QED radiation off the leptons. The 
package PHOTOS [34] 3 is used to handle multiple photon 
emission, with enforced p± -ordered radiation. 



q a i- q a j- 




1 1+1 i+ 



Fig. 2 Examples of mixed QCD0EW contributions included in the 
POWHEG BOX assuming factorization. 

In conclusion, let us note that the inclusion of all the the- 
oretical elements as described above allows to obtain pre- 
dictions for dilepton production in hadronic collisions tak- 
ing into account the interplay between QCD and EW con- 
tributions in a single computational framework. In particu- 
lar, the matching of the complete NLO SM corrections with 
the showers of both colored particles and photons enables to 
perform simulations of the differential cross sections in the 
presence of mixed QCD®EW effects, in addition to the con- 
tribution of higher-order QCD and QED corrections. In the 
absence of a complete calculation of 0(a em a s ) corrections 
to DY processes 4 , the tool enables to obtain, using factoriza- 
tion, leading logarithmic predictions for mixed QCD<g)EW 
contributions, as schematically shown in Fig. 2. 

3 Phenomenological results 

In this Section, we show the results of a variety of tests per- 
formed to validate the new theoretical and computational 

3 We use PHOTOS with leading log kernels, i.e. without process- 
dependent matrix element corrections already included in our approach 
in the exact NLO calculation. 

4 Partial results exist in the literature and are given by the one-loop EW 
corrections to W/Z+jet production at finite transverse momentum [35- 
38], NLO QCD corrections to the W/Z + y process [39^1] and two- 
loop virtual 0(cc em a s ) corrections to the NC DY [42]. However, the 
complete and non-trivial combination of all the above substructures is 
still unavailable. 



ingredients implemented in the POWHEG BOX. We first 
present comparisons against the benchmark predictions of 
the HORACE generator [11, 12] at NLO EW accuracy 5 . 
Then we address a phenomenological analysis of the in- 
terplay between QCD and EW corrections of interest for 
physics studies at the LHC. 

3.1 Input parameters and event selection 

We considered the lepton pair hadroproduction process pp — > 
Y,Z — > + (X), at the cm. energy y/s = 14 TeV. We 

used the MRST2004QED set of PDFs [43] with renormal- 
ization and (QCD and QED) factorization scale Hr = jj,f = 
M z for the comparison with HORACE and the QCD<K>EW 
simulations around the Z resonance, and )J,r = jip = M/+/- 
(the lepton pair invariant mass) for the numerical results on 
the combination of QCD and EW corrections well above the 
peak. For both the final-state leptons, we applied the follow- 
ing cuts on the transverse momentum and pseudorapidity: 

pf>20GeV, |%fc|<2.5, (2) 

which approximately model the acceptance of the ATLAS 
and CMS detectors at the LHC. In addition we also applied 
a cut on the invariant mass of the lepton pair of M/+;- > 
50 GeV. We considered, for simplicity, "bare" leptons, i.e. 
nearby photons are not recombined with the leptons. The re- 
sults have been obtained in the a e m(0) scheme, using the in- 
put parameters listed in Table l 6 , which coincide with those 
adopted in the tuned comparisons of EW predictions for Z 
boson observables described in Ref. [44]. The overall choice 
of input parameters, PDF set and acceptance cuts allowed 
us to check first that the NLO EW predictions of both HO- 
RACE and the new POWHEG version perfectly agree with 
the benchmark results for the integrated cross sections given 
in Ref. [44] . For the coupling of the external photons to the 
charged particles needed for the evaluation of the photonic 
corrections we used a eni = a em (0), and a s (Mz) — 0.1205 
with NLO evolution for the results about the combination of 
QCD and EW contributions. 

3.2 Comparisons with HORACE 

The results of the validation of the pure NLO EW correc- 
tions in the POWHEG BOX framework are shown in Figs. 3 
through 6 by comparisons with HORACE simulations. The 
invariant mass region around the Z resonance and well above 
it (M/+I- > 1 TeV) are considered in Fig. 3 - Fig. 4 and Fig. 5 

5 Note that HORACE implements completely independent EW form 
factors and real photon matrix elements, computed in the mass regu- 
larization scheme. 

6 The CKM mixing matrix is set to the identity matrix in the calculation 
of EW loop corrections. 
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Table 1 The input parameters used in the numerical simulations. 



Ofe m (0) = 1/137.03599911 
r z = 2.4924 GeV 
r w =2.141 GeV 
m e = 510.99892 KeV 
m„ = 0.06983 GeV 
m d = 0.06984 GeV 



M z = 91.1876 GeV 
sin 2 9w = l-M^/M z 
a s (M z ) = 0.1205 
m M = 105.658369 MeV 
m c = 1.2 GeV 
m s = 150 MeV 



M w = 80.37399 GeV 
M Higgs = 125 GeV 

m T = 1.77699 GeV 
m, = 174 GeV 
m A = 4.6 GeV 
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Fig. 3 Upper panel: lepton-pair invariant mass distribution according 
to POWHEG BOX and HORACE at NLO EW accuracy, in the reso- 
nance region. Lower panel: relative deviations, in percent, between the 
predictions of the two generators. 



Fig. 5 The same as Fig. 3 in the region M ; + ; - > 1 TeV. 
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Fig. 6 The same as Fig. 5 for the lepton transverse momentum distri- 
bution. 



Fig. 4 The same as Fig. 3 for the lepton transverse momentum distri- 
bution. 



- Fig. 6, respectively. For these comparisons we took the 
limit a s — > numerically in POWHEG. In the upper pan- 
els, we show the NLO predictions of the two codes for the 
invariant mass and p l ± distributions, together with the refer- 
ence Born results. The lower panels show the relative differ- 
ence (in percent) between the two programs, where the error 
bars correspond to la numerical uncertainties. 



We observe that the results of the two programs are in 
good agreement, both in the peak region and in the high 
tails. This means that all the EW ingredients have been cor- 
rectly included in POWHEG. In particular, the agreement 
in the Z resonance region is a test of the correct calculation 
of the NLO virtual and real contributions, as well as of the 
generalized subtraction procedure and treatment of collinear 
photon singularities. Indeed, the EW corrections in the peak 
region are known to be dominated by the mechanism of fi- 
nal state photon radiation, which receives contributions from 
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all the above components. On the other hand, the NLO EW 
corrections to the tails of the distributions give rise to large 
negative contributions due to Sudakov-like logarithms of the 
form a em \og(s / My) . The latter come from the exchange of 
weak gauge bosons in the loops, which therefore display in 
POWHEG the expected behavior at high energies 7 . 

In summary, the implementation in the POWHEG BOX 
of all the NLO EW corrections to the NC DY can be consid- 
ered fully under control. 



3.3 Combined effect of QCD and EW corrections 

The full results of the POWHEG BOX for the combined ef- 
fect of QCD and EW corrections to the invariant mass and 
lepton transverse momentum distribution are shown in Fig. 7 
- Fig. 8 in the resonance region, and in Fig. 9 for Mj+ t - > 
1 TeV. Further numerical results concerning the transverse 
momentum distribution of the Z boson in the resonance re- 
gion are shown in Fig. 10 - Fig. 13, in comparison with AT- 
LAS [49] and CMS [50] data at y/s = 7 TeV. 

For Fig. 7 - Fig. 9 the complete predictions have been 
obtained by matching the NLO QCD and EW corrections 
with QCD (PYTHIA version 6.4 [46]) and QED (PHOTOS) 
showers. For the sake of comparison, the pure QCD pre- 
dictions of the standard POWHEG BOX are also given, to- 
gether with the pure NLO EW results. The absolute predic- 
tions for the various distributions are shown in the upper 
panels of each plot. The lower panels display the relative 
difference, in percent, between the results of the new ver- 
sion of the POWHEG BOX and the standard QCD release, 
as well as the relative effect due to the genuine NLO EW 
corrections. Therefore the comparison between the two lines 
in each lower panel of Fig. 7 - Fig. 9 provides a measure of 
the QCDtgEW factorization and, more precisely, of mixed 
leading logarithmic corrections at the order a'^a" ,m,n> 1. 

Generally speaking, one can notice that both QCD and 
EW corrections (and their combination) are necessary for a 
proper control of the normalization and shape of the distri- 
butions. Particularly, the invariant mass distribution in the 
resonance region receives large corrections and a significa- 
tive shape modification from QED corrections, as empha- 
sized in previous studies [12, 13, 45]. As shown in Fig. 7, the 
left tail of the distribution is enhanced by the mechanism of 
final state photon radiation by several tens of percent, while 
the peak value is reduced by about 20%. For this distribution 
the effect of the EW corrections largely exceeds that of QCD 
radiation, and the impact of mixed QED<g)EW contributions 
is substantial, especially in the left tail of the distribution. 
The same kind of effect is present in the forward-backward 



7 Note that we do not take into account in our calculation the con- 
tribution due to real radiation of massive gauge bosons, which is 
known [47, 48] to partially compensate the negative virtual Sudakov- 
like corrections. 
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Fig. 7 Upper panel: lepton-pair invariant mass distribution around 
the resonance according to the full QCD&EW predictions of the 
POWHEG BOX (PWG EW + PYTHIA+PHOTOS), the standard QCD 
POWHEG BOX (PWG + PYTHIA), the LO and the NLO EW approx- 
imations. Lower panel: relative difference, in percent, between the full 
QCD®EW predictions and the pure QCD ones (green, dotted line), in 
comparison with the relative contribution due to pure NLO EW correc- 
tions (red, solid line). 
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Fig. 8 The same as Fig. 7 for the lepton transverse momentum distri- 
bution. 



asymmetry as a function of the lepton pair invariant mass 
for M/+[- below the Z mass, as we checked explicitly in our 
simulations. For the lepton p± the well-known overwhelm- 
ing QCD effects are by far dominant over the EW contri- 
butions, whose shape is washed out by QCD radiation (see 
Fig. 8). For such a distribution, the interplay of the particu- 
larly large QCD corrections with the ten percent level EW 
effects gives rise to mixed contributions of the order of sev- 
eral percents close to the peak, as clearly visible in the lower 
panel of Fig. 8. 

Non-negligible QCD(g)EW corrections are also present 
in the very high tail of the invariant mass distribution shown 
in Fig. 9. In this region the large EW corrections, enhanced 
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Fig. 9 The same as Fig. 7 in the region M/+ ( - > 1 TeV. 
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Fig. 10 The normalized differential cross section as a function of p\ 
for the full range up to 350 GeV (whole plot) and in the range //[ 
< 20 GeV (inset). The ATLAS data (for bare muons) are compared 
with the predictions of the POWHEG BOX with full QCDigiEW con- 
tributions and QCD corrections only, and according to two PYTHIA 
versions. 



Fig. 11 Upper panel: ratio of the POWHEG BOX QCD predic- 
tions over the ATLAS data, according to two different version of 
PYTHIA. Lower panel: relative difference, in percent, between the full 
QCDigEW predictions and the pure QCD ones of POWHEG BOX for 
two different PYTHIA versions. 
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Fig. 12 The same as Fig. 10 for CMS data (for bare muons) up to 
600 GeV. 



by Sudakov-like logarithms, in association with QCD radi- 
ation induce mixed contributions which grow from a few to 
several percents, as can be seen in the lower panel of Fig. 9. 

Last but not least, we show for illustrative purposes in 
Fig. 10 - Fig. 13 our results for the Z transverse momen- 
tum distribution, whose knowledge is crucial for a precise 
measurement of the W mass and whose description repre- 
sents a challenge for QCD. The predictions of the POWHEG 
BOX are compared with ATLAS data [49] (see Fig. 10 - 
Fig. 11) and CMS data [50] (see Fig. 12 - Fig. 13) at y/s = 
7 TeV 8 . For these results, we imposed the acceptance and 
invariant mass cuts quoted in the experimental papers, we 
focused on the data referring to bare muons and we used 



We used the data available at the HepData repository 
http://hepdata.cedar.ac.uk. 



the MSTW2008 NLO PDF set [52]. Because of the con- 
siderable dependence of the /r[ distribution on adjustable 
non-perturbative parameters, we considered POWHEG in- 
terfaced to two different versions of PYTHIA generator, na- 
mely PYTHIA version 6.4 [46] and version 8.1 [51]. For 
each PYTHIA version, we used the default values for the pa- 
rameters describing the underlying event and shower effects. 
While higher-order QED radiation is treated with PHOTOS 
when POWHEG is interfaced to PYTHIA 6.4, the simula- 
tion performed with PYTHIA 8.1 uses its internal available 
tool for QED shower. 

In Fig. 10 and Fig. 12 we show the absolute predictions 
of POWHEG for the normalized differential cross section 
as a function of p\ in comparison with ATLAS and CMS 
data, respectively. The POWHEG results are presented both 
with and without NLO EW and QED Shower corrections, 
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Fig. 13 The same as Fig. 1 1 for CMS data. 



and according to the two aforementioned PYTHIA versions. 
The inset figures show the data-theory comparison in the 
low region. Albeit an optimal description of the Z trans- 
verse momentum distribution over the full range seems to be 
problematic, there is a rather good agreement with the data, 
especially when considering the POWHEG BOX interfaced 
to PYTHIA version 8.1. The latter statement can be under- 
stood by looking at the upper panels of Fig. 11 and Fig. 13. 
They show the ratio of the POWHEG BOX QCD predictions 
over the ATLAS and CMS data, as obtained with PYTHIA 
version 6.4 and version 8.1, respectively. Especially in the 
low p\ region, POWHEG interfaced to PYTHIA 8.1 pro- 
vides a much better description of the data, both for ATLAS 
and CMS. The lower panels of Fig. 1 1 and Fig. 13 show in 
conclusion the impact on the Z transverse momentum dis- 
tribution due to NLO EW and multiple photon corrections, 
as well as to 0(a em a s ) contributions. As can be noticed, 
these contributions are largely independent of the particular 
version of the QCD/QED Shower generator under consider- 
ation. They introduce a negative correction of some percents 
in the first bins of the jr\ distribution, while yielding a mod- 
erate effect of about +2% for p\ around 20 - 30 GeV. 

As a concluding remark, it is worth noting the the pos- 
sibility of simulating the contribution of mixed QCD<g)EW 
corrections with the new tool is due to the particular fac- 
torized form of the EW corrections in the QCD POWHEG 
framework. In this respect, our predictions contain comple- 
mentary information with respect to the calculation of Refs. 
[17-20], where QCD and EW corrections are combined ad- 
ditively. Moreover, the approach includes details of the in- 
terplay between QCD and EW radiation which are not taken 
into account in the experimental procedure typically followed 
at the LHC of applying EW ^-factors to the predictions of a 
reference QCD program. 



4 Conclusions 

We have generalized previous work on the inclusion of the 
EW corrections to single W production in the POWHEG 
BOX to cope with the NC DY process. We have added NLO 
EW and QED multiple photon corrections to the native NLO 
and PS QCD contributions. To this aim, we have exploited 
the process-independent structure of the POWHEG BOX 
framework and resorted to various general techniques al- 
ready successfully developed for W -boson production. 

We have provided evidence of the accuracy of the ap- 
proach and presented a sample of phenomenological results 
about the combination of QCD and EW corrections to lep- 
ton pair production in hadronic collisions at LHC energies. 
In particular, we have shown that the leading contribution of 
mixed QCDtgEW corrections may affect some observables 
at the level of several percents, beyond the separate effect 
of strong, weak and electromagnetic corrections. Therefore, 
the exact calculation of 0(a em a s ) corrections would be de- 
sirable for a better theoretical control of the DY cross sec- 
tions. 

The new tool enables to obtain precise and realistic pre- 
dictions for the NC DY process, in the same way as the 
charged-current channel. It is available at the web site of 
POWHEG BOX. Further theoretical ingredients, such as pho- 
ton-induced processes [12, 19, 53, 54] and presently ne- 
glected higher-order corrections, will be made available in 
future releases. 

We expect that the results here presented will facilitate 
the work of experimental analysis and data interpretation at 
hadron colliders and will allow in the future a rather simple 
inclusion in the POWHEG BOX of EW effects to further 
processes of physics interest. 
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